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Abstract 

Recently, a class of gravitational backgrounds in 3 + 1 dimensions have 
been proposed as holographic duals to a Lifshitz theory describing critical 
phenomena in 2 + 1 dimensions with critical exponent z > 1. We continue our 
earlier work [17], exploring the thermodynamic properties of the "black brane" 
solutions with horizon topology M 2 . We find that the black branes satisfy the 
relation 6 = -^-fTs where £ is the energy density, T is the temperature, and 
s is the entropy density This matches the expected behavior for a 2 + 1 
dimensional theory with a scaling symmetry (xi, X2) — > X(x±, X2), t — > X z t. 



1 Introduction 



Since the Maldacena conjecture [1], holography has offered an interesting new tool to 
explore strongly coupled field theories (for a review, see [2]). In this framework, black 
hole backgrounds are dual to strongly coupled plasmas, and using these backgrounds, 
one can extract hydrodynamic and thermodynamic properties of the plasma [3-5]. 

Recently, much effort has gone into describing quantum critical behavior in con- 
densed matter systems using holographic techniques (for a review, see [6]). Quantum 
critical systems exhibit a scaling symmetry 

t — > X z t, Xi — > Xxi (1) 

similar to the scaling invariance of pure AdS (z = 1) in the Poincare patch. From a 
holographic standpoint, this suggests the form of the spacetime metric 

ds 2 = L 2 (r 2z dt 2 + r^dx'dxHij + ^] , (2) 



r 



where the above scaling is realized as an isometry of the metric along with r — > A -1 r 
(for our purposes, i = 1,2). Other metrics exist with the above scaling symmetry, 
but also with an added Galilean boost symmetry [7-9] which we will not consider 
here (black brane solutions in these backgrounds were discussed in [10]). There has 
also been some success at embedding a related metric into string theory [11] with 
anisotropic (in space) scale invariance. This may serve as a template for embedding 
metrics of the form ([2]) into string theory. Here, however, we will continue our study 
of the model in [12] (some analysis of generalizations of this model appear in [13-15]). 
The authors of [12] constructed a 4D action that admits the metric ([2]) as a solutions, 
which is equivalent to the action [13] 

S = — ^- f d'x^-g (r-2A- \T, W T^ - ^A^ M ) (3) 
[T = dA) with terms in the action parameterized by 



L 2 L 2 

The solution discussed in [12] has the metric ([2]) and 



2Z A= _lg+g+4 (4) 



2z(z- 1) 



A = L'-^^j^dt (5) 



lr This background was earlier studied in [16], however, without an action principle that admits 
the above metric as a solution. 
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with the identification z = Z and L = L or the identification z = 4/Z, L = and 
is defined for z > 1 (with z = 1 giving AdS^. 

In our current work, we will analyze the thermodynamics of black brane solutions 
which asymptote to (j2J) 0- These black brane solutions have been studied numerically 
in [17-19]. However, we will show that their energy density £, entropy density s and 
temperature T are related by 

£ = -^—Ts (6) 

2 + z v ; 

using purely analytic methods. 

In the following section we will introduce the ingredients needed to prove (jHJ). 
Particularly important is the existence of a conserved quantity, used to relate horizon 
data to boundary data. In the final section, we combine these ingredients into the 
result (jSJ). We also show that this relation is expected in the dual field theory as a 
result of the scaling symmetry (xi,x 2 ) ~" * A(xi,X2), t — > \ z t. 

2 Summary of earlier results 

In our earlier work [17], the action (jHD was reduced on the Ansatz 

ds * = ^Mr) dt 2 + e 2B(r)^ dx l ) 2 +{dx 2 ) 2 )+e 2C(r) dr 2 

A = e G ^dt (7) 
to give a one dimensional Lagrangian 

L 1D = Ae^ 2B+A -^dBdA + 2e^ B+A - c \dB) 2 + ^- A+2B - c+2G \dG) 2 

-2Ae {A+2B+c) + - c 2 e (~ A + 2B + c + 2G ), (8) 

The equations of motion following from this action have solutions given by (J2J, (151). 
Further, there are the known black brane solutions that asymptote to AdS4, 



A J V r 2 f(r) 



f(r) = 1-5 



r 



3 



^ = 0. (9) 



2 We constructed numeric solutions for this system in [17] (seen as the large tq limit of this earlier 
work), however, our current work does not depend on any numeric analysis. 
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For the remainder of the paper, we will be concerned with black branes which asymp- 
tote to the Lifshitz background (j2J), (EJ) for z ^ 1. 

In [17], a Noether charge was found which is associated with the shift 



(A{r)\ 
B(r) 
C(r) 

\G(r)J 



fA(r) + 5\ 
B{r) - I 
C{r) + 

\G{r)+5j 



(10) 



with 5 a constant. The above represents a diffeomorphism which preserves the volume 
element dtdx l dx 2 . This is why it is inherited as a Noether symmetry in the reduced 
Lagrangian. The associated conserved quantity is 



( 2e (A+2B-C) 9A _ 2e{A+ 2B-C) dB _ e{ - A+ 2B-C + 2G) dG) = ^ 



'111 



2.1 The perturbed solution near the horizon 

We begin by first reviewing the results found in [17] for the expansion near the 
horizon. We require that e 2A goes to zero linearly, e 2C has a simple pole, and e G 
goes to zero linearly to make the flux dA go to a constant (in a local frame or not). 
Further, we take the gauge B(r) = ln(Lr) for this section. We expand 

A(r) = In (r z L (a$(r — r )z + a ai(r — r ) i + ■ • • j j , B(r) = ln(rL) 
C(r) = In (- (co(r - r )"5 + a{r - r )2 + •■•)], (12) 



, . / L 2 r z 2z(z-l) , , . , , 2 

G(r) = In I —d — (a g {r - r ) + a gi{r - r ) +■■ 

Note that by scaling time we can adjust the constant by an overall multiplicative 
factor (note the use of in the expansion of G(r) as well, as e G multiplies dt for 
the one- form A). We will need to use this to fix the asymptotic value of A(r) to be 
exactly \n(r z L) with no multiplicative factor inside the log. 

We plug this expansion into the equations of motion arising from (jSJ), and solve 
for the various coefficients. We find a constraint on the 0th order constants: as 
expected not all boundary conditions are allowed. We solve for c in terms of the 
other g and r , and find 

e _ y/(2z + g*r (z-l))r$ 
° ^(z 2 + z + 4)r 2 
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All further coefficients (e.g. g\) are determined from the two constants tq and go- 
We evaluate (TTTjl at r = r using the above expansion to find 

D = r ° +3jL2a ° . (14) 
c 

This must be preserved along the flow in r. We will use this to relate constants at 
the horizon to coefficients that appear in the expansion at r = oo. 

2.2 The perturbed solution near r = oo 

We now turn to the question of the deformation space around the solution given in 
(GO) and (OS). We take the expansion of the functions 

A(r) = hx{r z L) + eA x {r), B(r) = ln(rL) + eB^r) 
C(r) = In f^j +eCi(r), G(r) = In r * 

Using straightforward perturbation theory, we may find the solutions in the Bi 
gauge 

= Co ( ^ ~ 2) r~ z ~ 2 + C 2 (2 2 + 3^ + 2 - (z + 1) 7 ) r-f-^i (16) 

I 2 + 

5i(r) = (17) 
d(r) = -C (z-l)r- 2 - 2 + C 2 (z 2 -7z + 6-(^-l) 7 )r-i- 1 -i (18) 

d(r) = C o 2( f + + 2 2) r """ 2 + C24z(z + ljr-i-^i (19) 
where we have defined the useful constant 




7 = - 20^ + 20. (20) 

In the above, we have dropped certain terms in the expansion from [17]. We have 
dropped them so that we meet the criterion of [20] to be sufficiently close to the 
Lifshitz background (j2J), ©. We may evaluate the conserved quantity, and we find 

D = _2( S -l)( S -2)( S + 2)£,' 

z 

One may worry that nonlinearities may contribute to the value of this constant. 
However, one may examine the powers of r available in [17], and quickly be convinced 
that the higher nonlinear contributions will be zero once we meet the criterion of [20] 
(i.e. dropping the "bad" modes). Hence, the above constant is the value of Do 
throughout the flow in r. 
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2.3 Gauge invar iance 



In the previous sections, we have gauge fixed by taking -Bi(r) = 0. Here, we write 
down the linearized gauge transformations that will allow us to switch to other gauges 
in perturbation theory (used near r = oo). The transformation 

A 1 (r) -+Ai(r) + -5(r), B^r) B 1 (r) + -5(r) (22) 

Cx(r) - d(r) + rd r ( S -^-) , G 1 (r) -> G x (r) + -5(r) 

\ r J r 

corresponds to infinitesimal coordinate transformations r — > r + e<5(r). One can see 
that such a shift leaves the equations invariant to leading order in e when expanding 
about the solution (j2D, (JSD- 



3 The black brane thermodynamics 

In the above, we have calculated the conserved quantity D in two regions: near the 
horizon r = r and in the asymptotic region r = oo. We may use this to solve for C 
in terms of the horizon data 

= I r z +3 a z 
2co(z-l)(z-2)(z + 2)' 1 1 

We will see that Co is proportional to the energy density 8 of the background. 

Indeed, the authors of [20] identified the energy density of the background in 
terms of the coefficient of the r~ z ~ 2 term in the expansion at infinity for any z. This 
mode was also identified in [17] as the mass mode using background subtraction. 
However, for 1 < z < 2 there were additional divergences that were not canceled. 
These were cured using the local counter terms of [20]. 

One may not, however, directly use the results of [20] because of a different choice 
of gauge: above we use £>i(r) = and the authors of [20] use C\(r) = 0. We may 
easily switch to this gauge by taking 

5( r ) = -r I ^^dr (24) 



and transforming the other fields appropriately (near r = oo). In the above integral, 
we make sure to take the constant of integration so that at large r the correction 
fields vanish. We find that in the C*i(r) = gauge 
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where ■ ■ • are the other terms we are not concerned about. This allows us to compare 
our results directly to the calculation of [20]. We identify 2Ai(r) = f(r) where f(r) 
appears in equation (5.27) [20]. Therefore, we identify 

ci,rs = -(z-1)C (26) 

where c^rs appears in [20] as the coefficient of r~ z ~ 2 . Therefore, equation (5.31) 
of [20] becomes 

4(s-2)(z-l) „ 2rg +3 a 

8 ~ z C ° ~ j2 + zfo' (27) 

where £ is the energy density. However, we should note that this is a unitless energy. 
Restoring units we find 

= 2rg+ 3 a 1 

(2 + z)c 1QttG 4 L' y ' 

From the metric, it is easy to read off the area of the horizon, and therefore the 
entropy density 

107rG4 

Further, one can easily read off the temperature from the expansion at the horizon 
[17] 

- r o +1 «o (30) 



From this, we may read an interesting thermodynamic relationship 



£ = ih Ts ' (31) 



which is the main result of our work. 

We compare this expression with the known z = 1 black brane solution in AdS^. 
For this we have £ = 167T q 4L , T = and s = 167T q 4 ■ These satisfy the relations 
d£ = Tds,£ = \Ts. The second relation agrees with our expression above for z — 1. 

In general, one may use 2 relations of the form 

d£ = Tds (32) 
£ = KTs (33) 

(K a constant, which for our purposes is a function of z) to find the functional forms 
£(r ), £{T) and s(r ), s(T). First, one may use (1331) to eliminate ds in the relation 



6 



(1321) . directly relating dE and dT. One may integrate this to find S(T). One may 
then use this in ( 1331) to find s(T) = ■ Such a procedure will furnish T(r ) and 
so we can find S(r ) and s(r ). Doing so, we find 



l-A - 
K 



^ = 16^? fW = e (ra) ' r(r »»=(l6^J (34) 

where O is a constant of integration independent of tq. has units of length i--^ so 
that the units of s and 8 are canonical. We can rewrite in the following way. tq is 
a metric parameter, and the temperature should only depend on metric parameters, 
not IGttG^. Further, the energy density should depend on WttG^, with one power of 
this factor in the denominator. Therefore, we can say that O oc 16 J G4 ■ To make up 
the rest of the units, there is only one dimensionful parameter that one can use: L. 

K 

Therefore, we take that O = 1fi ~ r n(z). This gives 



K(z) 



l-K(z) 

1 f K(z)4irrl\ *W 
n(z) 



where we now explicitly write that K is a function of z \K(z) = ^tj) • To the 
right of each of these expressions, we have given the relation in terms of only the 
thermodynamic variable T, rather than referring to the geometric variable r Q . 

One may use the above relations to show that the pressure P is related to the 
energy density by 2P = zS, which matches the holographic result of [20]. Actually, 
starting from the equation of state 2P = z£ one can derive the relation fl3TT) . 

The value of n(z), which is a unitless number, may be determined numerically 
by plotting log(LT) vs log(r ), but we did not attempt to do this here. 

The slope of the graph of log(T(ro)) vs log(ro) is 

log(T(r )) = — 7-^— log(r ) + constant = zlog(r ) + constant, (36) 

K[z) 

which can be quantitatively checked using our earlier data [17]. 

We may easily compare the thermodynamic relationship (j3TT) with the expected 
behavior from a system with a scaling symmetry (x\, 22) — * A(xi, X2), t — > X z t. This 
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scaling symmetry (along with SO (2) symmetry rotating the Xj) implies that the 
dispersion relation is 

u 2 = a\k\ + klY = a 2 k 2z (37) 

where a is some parameter to restore canonical units. We will work with a finite 
system (a box with sides of length £), and assume that the occupation number of a 
given mode (with energy E n ) in the box is Q' (e _/3Sn ) e~^ En . This way, the energy 
of the system is written as 

Inside the box, hi = jrTii, and we have d = 2 spatial dimensions, so there are two 
rij. We approximate the sum by an integral, and realize that the density of integer 
(ni,n 2 ) lattice points is uniform in the (ni,n 2 ) plane to find 



Q = 2tt J n l --Q [eT^K- n ) j , n { + n\ ee n 2 . (39) 
Redefining the integration variable, we find 

02 poo 

- 1 dXX^Q (e" A ) ee Z -~V^, (40) 

where the constant H is independent of /3, and we have assumed that the function Q is 
well behaved at infinity (this is used to construct S). This yields all thermodynamic 
quantities 

and so indeed relation (!3T|) holds, as well as all subsequent formulae. Using the above 
argument, we can generalize the result to arbitrary spatial dimension d where the 
relation ( 1311) becomes 

8 = Ts-^—. (42) 



d + z 



which generalizes the function K(z) K(z,d) = (we also promote n(z 



d+z 

n(z, d)). For d = 3 and z = 1 (giving K — |) one can check that fT4"2"|) and flHol) are 
correct by comparing to the known results for black D3 branes [3] for the functions 
s(T) and S(T) up to the normalization n(z,d). The functions of s(r ), S(ro), T(r ) 
need to be modified by promoting Tq — > Tq, so that the "area" is measured appropri- 
ately. With this modification, the expressions for s(r ), £(r ), T(r ) also agree; see 
for example [21], where similar coordinates are used. 
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